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Abstract 

A system of nonlinear integral equations equivalent to the steady- 
state Krook kinetic equation was used to model the flow from a low-thrust 
axisymmetric nozzle. The mathematical model was used to numerically 
calculate the number density, temperature, and velocity of a simple gas as 
it expands into a near vacuum. With these quantities the gas pressure and 
flow directions of the gas near the exit plane were calculated and compared 
with experimental values for a low-thrust nozzle of the same geometry and 
mass flow rate. 


m Introduction 

00 

r» 

£ A coordinated program to study gas flows from low-thrust nozzles 

o (*1N) is in progress at NASA-Lewis Research Center. 1 This effort includes 

numerical mathematical modeling and experimental measurements to 
° determine plume parameters. 

\ In two previous papers, model results were presented for plume 

2 flows from a two-dimensionaP nozzle and a three-dimensional 

axisymmetric 3 nozzle. No comparison of the model results were made with 
experimental data. In this paper, for the first time, the model results and 
experimental data are compared. The two earlier papers showed the 
feasibility of utilizing the Boltzmann 4 equation with a Krook-type relaxation 
5 scattering term for binary collisions to model plume flows for low-thrust 
nozzles when the flow goes through the transition region (continuum to free- 
molecular flow). In particular, it was demonstrated that an iterative process 
for solving the Boltzmann equation would converge to give a steady-state 
solution using a reasonable amount of computer time. The spherical 
computational volume was quite small (-3 cm radius) and the nozzles at the 
exit plane were only 3/4 cm in diameter. 

The model work described here is for a larger nozzle (exit plane 
diameter -3.2 cm) and extends the computational volume (-20 cm radius) 
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in an attempt to predict parameters in the near and intermediate regions of 
the plume. The nozzle geometry and propellant flow rate in the 
experimental setup are identical to those assumed in the model, thus model 
predictions can be directly compared with the experimental measurements. 

Kinetic Theory Model 


Mathematical Formulation 

The Krook kinetic equation 5 or iterate equation 4 for steady-state 
flow that was solved numerically by integrating over velocity space along 
a characteristic line is 

v V,f(v,/) =A(nW.-f) (1) 

where t, (shown below) is an equilibrium Maxwellian velocity distribution 
normalized to one at a given r. 

f. = {1/(2Ttk7) 3 *}exp-{(m/(2k7)(v'-u) *} (2) 


The boundary conditions for f (probability distribution function) are 
discussed in the next section of this paper. 

The first iteration values for n(r), u(i), and T{r) (number density, 
mean velocity, and temperature) at each grid point can then be calculated 
using the following definitions 


n(r)=i II f(r,v)d 3 v 

(3) 

u{t) = (1/ n) JJ / irf (r,v) d 3 ir 

(4) 

T[i) = (m/(3k Tn)} J J J (v-i/) 2 f(r,v) tfv 

(5) 


The values of n, u, and T are put into Eqs. (1) and (2) and the 
process 3 is repeated until convergence is obtained. At convergence the 
values for n, u, and T at the grid points do not significantly change from 
the values obtained in the previous iteration. Because the nozzle is 
axisymmetric, only parameter values in one plane through the axis of 
symmetry are necessary. Figure 1 shows the plane polar grid system used 
in the model calculations. 

Nozzle Parameters and Flowfield Boundary Conditions 

As mentioned earlier the experimental flow conditions and nozzle 
geometry were identical with those in the model. The nozzle was 
axisymmetric with a half angle at the exit plane of 20 deg. The nozzle throat 
diameter was 3.18 mm and the area ratio (exit/throat) was 100. The plume 
gas was N 2 with a flow rate of 6.8 x10 s kg/s. 

To simulate the 0.02 Pa experimental background pressure in the 
vacuum chamber, the plume in the model was assumed to expand into a 
background number density set at 1.2x10'* N 2 molecules/m 3 with a 
temperature of 300 K at zero mean velocity. This was accomplished in the 
model by assuming experimental vacuum chamber conditions for the 
surface located at 200 mm with angles from the center line greater than 50 
deg (See Fig. 1). For the surface at 200 mm in the first 50 deg from the 
centerline, the flow was assumed to have only an outward velocity 
component. 

The flow parameters at the nozzle exit plane used for the inflow 
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boundary conditions were obtained from Boyd 8 - 7 who utilized a discrete 
simulation Monte Carlo (DSMC) computational method to model the flow 
inside the nozzle. Boyd’s model assumes diffuse scattering along the inside 
of the nozzle wall. Thus the boundary layer effect on the flow external to the 
nozzle will be inherent through the exit surface boundary conditions. The 
parameters at exit surface grid points are given in Table 1. Boundary 
conditions for the external wall of the nozzle were modeled with f =0. 

Convergence Criterion and Computer Time 

The convergence to a solution for the iteration process was 
assumed when changes in the quantities n, u, and T compared to previous 
iteration values approached zero. A solution for the 539 grid points shown 
in Fig. 1 took approximately 20 h of CPU time on the CRAY -XMP at the 
NASA-Lewis Research Center. An equivalent solution was obtained with 
about 1/2 the number of grid points and correspondingly less computer 
time. Most of the computer time was taken in the relatively high-number 
density regions near the nozzle exit plane and along the radii lines at large 
angles (angle > 50 deg) where the number density is near the background 
level. 

For convergence at a grid point, ctest, defined to be the square 
root of the sum of the squares of the fractional changes in the number 
density, the components of velocity, and the temperature was required to 
be small. Typically the convergence criterion for ctest was <0.01. A 
correction factor that ranged from 1 .0 to 1.1 as r in Fig. 1 ranged from 46.4 
mm to 200 mm was applied to the number density to insure that the mass 
flow rate at each radius was maintained at - 6.8x1 0 -5 kg/s. This correction 
was needed when the grid point separations were large. 

Experimental Data 

The experimental data used in this paper were taken in a vacuum 
tank for space-simulated tests at NASA-Lewis Research Center . The 
details of the experimental setup are described in a paper by Penko 8 et al. 
and the correction factor that was applied to the pressure measurements 
due to the shock wave caused by the measuring probe is discussed by 
Shapiro.* A conical probe was used to measure the flow angles and Pitot 
tubes with diameters of 1 mm and 6.4 mm were used to measure the total 
pressure in appropriate regions of the plume. The 1-mm tube was used 
near the exit plane where large pressure gradients exist and the 6.4-mm 
tube was used in the far-field region of the plume. An overall measurement 
accuracy in the pressure measurements of ±9% was estimated at the exit 
plane. This improved to ±2% in the far field. The 30-deg conical probe 
used to measure the flow angle in the farfield was 6.4 mm in diameter at the 
base and had an estimated accuracy of ±1 deg. 

Model Results Compared with Experimental Data 

Model Prediction Overview and Parameter Manipulations 

The model results are depicted in contour plots for the grid points 
as seen in Fig. 1 with the exclusion of the grid points that have r values 
less than 46.4 mm. Figures 2-5 are contour plots for the number density, 
axial component of velocity, radial component of velocity, and temperature, 
respectively. These plots give a complete visual perspective of the model 
results. The contour plot in Fig. 6 shows the Pitot pressures as calculated 
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from the model which can be compared directly with the data. 

Pitot pressures were calculated from the model values of n, u, 
and T at the grid points with shock wave and supersonic rarefied flow 
effects built into the model results. Since the Pitot probe acts as an 
intrusive device around which a shock wave develops in the supersonic 
flow, the pressure drop across the shock wave must be calculated in order 
to compare with the experimental data. Thus, the Rayleigh supersonic Pitot- 
tube formula 9 was applied to the model prediction for pressure. 

A second modification to the model prediction for pressure was 
necessary due to the rarefied nature of the flow. The following 
equation'- 10 was applied to the model results 

log( P OT /PJ = .089 -. 1 20 log(Pe p ) for Re„<, 5.6 (6) 

P m is the corrected model predicted Pitot pressure which is to be 
compared directly with the experimental data. Po, and Pet, are, respectively, 
the model calculated values for the pressure behind the shock wave and 
the Reynolds number. The Reynolds number, Pq, , was calculated utilizing 
the following relation 


Pe p = p.U.D P 4i (7) 

where p. and U. L are the gas density and velocity of the freestream 
preceding the shock, whereas D„ and |i are, respectively, the probe 
diameter and the viscosity at stream temperature, T„ behind the shock 
wave. The temperature, Ty, was computed from the normal-shock relation 
for static temperature (See Ref. 9 , p. 1 18). Once, T„ was known, then the 
viscosity, p, was determined from the power-law relating viscosity of gases 
(See Ref. 1 , p. 6) by the following formula 


P = P™.(VU" (8) 

where p* is the value of the viscosity at the reference temperature, T ral , 
and n = 0.75. 


Direct Comparison of Flow Angle and Pressure with Experimental 
Data 

To compare the model predictions with the experimental data, the 
pressures and flow angles at the grid points in Fig. 1 were linearly 
interpolated to pressures and flow angles along the three radial lines at 
axial positions of 12 mm, 36 mm, and 120 mm as shown in Fig. 1 . 

Results of the comparisons at these three axial positions for 
pressure as a function of radial distance is shown in Fig. 7. Results for the 
flow angles are given in Fig. 8 for the 12-mm and 120-mm axial positions 
since there was no reliable experimental data at the 36-mm axial position. 
In general there is excellent agreement between the model predictions and 
the experimental data. The biggest variation between the model predictions 
and the experimental data is at the 12-mm axial position for both the 
pressure and the flow angle where the experimental data are not as 
precise. 
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Concluding Remarks 


A kinetic theory model was used to calculate flow parameters 
external to a low-thrust axisymmetric nozzle. The model predictions for flow 
angle and pressure were compared for the first time with experimental data 
taken for identical flow conditions. Agreement between the model 
predictions and the experimental data was very good, which gives one 
confidence that the model and numerical techniques used to calculate the 
flowfield parameters may be useful in calculating flows from other low- 
thrust nozzles. In particular, the model could be used to calculate the 
flowfield parameters under space conditions that would be prevalent around 
the space station being built by the United States. Preliminary model 
calculations using space conditions look promising. 
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Fig. 1 Grid system. Grid points are at intersections of the radii and circular 
lines. All grid points are in the y ■ 0 plane. The lines marked as 12 mm, 36 mm, 
and 1 20 mm are the radial lines along which experimental data was compared 
to the predictions of the model. The nozzle lip is at radius 46.4 mm and angle 
of 20°. 


Number Density Contour Plot 
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Fig. 2 Number density contour plot for flowfield shown in Flg.1. 
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Axial Valoctty Contour Plot 



Fig. 3 Axial velocity contour plot for 
flowfleld shown In Fig. 1. 


Tamperatura Contour Plot 
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Fig. 5 Temperature contour plot for 
flowfleld shown in Fig. 1. 


Radial Valoclty Contour Plot 
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Fig. 4 Radial velocity contour plot 
for flowfleld shown in Fig. 1 . 
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Fig. 6 Pressure contour plot for 
flowfleld shown In Fig. 1. 
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Comparision of Pitot Pressure 
with Model Prediction 




Prediction at 12 mm 


Prediction at 36 mm 


Prediction at 1 20 mm 

• 

Data at 12 mm 

• 

Data at 36 mm 


Data at 120 mm 


20 30 40 

RADIAL DISTANCE (mm) 


50 


60 


Fig. 7 Comparison of experimental Pitot pressures and model prediction at 
three axial positions. The last radial data point at each axial position 
corresponds to the end of the appropriate line In Fig. 1. 


Comparison of Flow Angle 
with Model Prediction 



Fig. 8 Comparison of experimental flow angles and model predictions at two 
axial positions. The last radial point at each axial position corresponds to the 
end of the appropriate line In Fig. 1. 
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TABLE 1 BOUNDARY CONDITIONS AT EXIT SURFACE GRID POINTS 


Angle 

Radial 

Axial 

Tamparatura 

Numbar 

Pitot 


Velocity 

Valoctty 


Density 

Pressure 

(Degrees) 

(nVa) 

(nVa) 

<K) 

(Molecules/ m 3 ) 

(Pa) 

0.0 

0.0 

1180 

99 

4.46 E +21 

272.0 

2.3 

35.0 

1170 

99 

4.51 E +21 

269.0 

4.6 

66.8 

1140 

102 

4.62 E +21 

262.0 

6.8 

89.6 

1120 

115 

4.30 E +21 

237.0 

6.8 

140.0 

1070 

148 

3.36 E +21 

173.0 

10.6 

174.0 

887 

201 

2.42 E +21 

116.0 

13.0 

215.0 

858 

278 

1.53 E +21 

61.2 

15.0 

233.0 

732 

338 

1.13 E +21 

36.9 

17.8 

237.0 

525 

421 

7.43 E +20 

16.4 

18.4 

216.0 

415 

467 

6.30 E +20 

11.0 

20.0 

135.0 

338 

536 

6.01 E +20 

9.2 






